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Abstract 

We present a measurement of the B° — > Apir~ branching fraction performed using the BABAR 
detector at the PEP-II asymmetric energy e + e~ collider. Based on a 232 million BB pairs data 
sample we measure: B(B° — ► Ap-K~) = [3.30 ±0.53 (stat.) ±0.31 (syst.)] x 10~ 6 . A measurement 
of the differential spectrum as a function of the di-baryon invariant mass m(Ap) is also presented; 
this shows a near-threshold enhancement similar to that observed in other baryonic B decays. 
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1 INTRODUCTION 



Observations of charmless three-body baryonic B decays have been reported recently by both the 
BABAR and Belle collaborations ^ ^ |3] ■ A common feature of these decay modes is the peaking of 
the baryon-antibaryon mass spectrum toward threshold. This feature has stimulated considerable 
interest among theorists as a key element in the explanation of the unexpectedly high branching 
fractions for these decays [U E] • We report a measurement of the branching fraction for B° decay 
to the Apir~ final state 5 . In the Standard Model this decay proceeds through tree level b — > u 
and penguin b — > s amplitudes. It is of interest to study the structure of the decay amplitude 
in the Dalitz plane and to test the afore-mentioned theoretical expectations. This channel may 
also be used to search for direct CP violation, and with the A hyperon in the final state, its spin 
self-analyzing weak decay to p ir, may be used, with increased statistics, to study the chirality 
structure of weak b — > s transitions |Hj and to probe T violation [UEj. 

2 THE BABAR DETECTOR AND DATASET 

The data sample consists of 232 million BB pairs corresponding to an integrated luminosity of 
210.3 fb~\ collected at the T(4S*) resonance with the BABAR detector. The detector is described in 
detail elsewhere [5]. Charged particle trajectories are measured in a tracking system consisting of 
a five-layer double-sided silicon vertex tracker (SVT) and a 40-layer central drift chamber (DCH), 
both operating in a 1.5-T axial magnetic field. A ring-imaging Cherenkov detector (DIRC) is used 
for charged-particle identification. A CsI(Tl) electromagnetic calorimeter (EMC) is used to detect 
and identify photons and electrons, while muons are identified in the instrumented flux return of 
the magnet (IFR). A BABAR detector Monte Carlo simulation based on GEANT4 jH] is used to 
optimize selection criteria and determine selection efficiencies. 

3 EVENT SELECTION 

We reconstruct A candidates in the A — > p?r decay mode as combinations of oppositely-charged 
tracks, assigned the proton and pion mass hypotheses, and fit to a common vertex. A fit to 
the invariant mass distribution of reconstructed candidates with a triple Gaussian gives RMS 
widths of 0.7MeV/c 2 , 1.3MeV/c 2 and 4.0MeV/c 2 for the narrow, intermediate and wide Gaussians 
respectively, with an average value of 1.8 MeV/ c 2 . Combinations with an invariant mass in the range 
1.111 - 1.121 GeV/c 2 are refit with a mass constraint to the nominal A mass |l()j . and combined 
with two additional tracks with opposite charges, each with momentum transverse to the beam 
greater than lOOMeV/c. 

Measurements of the average energy loss (dE/dx) in the tracking devices, angle of the Cherenkov 
cone in the DIRC, and energy releases in the EMC and IFR are combined to give a likelihood 
estimator for a track to be consistent with a given particle hypothesis. Tracks with likelihood 
ratios satisfying the very loose particle identification (PID) criterion L p /Lx > 1/3 and Lp/L^ > 1 
are assumed to be protons. In addition, the pion that originates from the B decay vertex must 
satisfy a loose PID criterion L n /Lx > 0.22 and L n /L p > 0.02. A Kalman fit to the full decay 
sequence is used to reconstruct the B vertex; only candidates with a fit probability P vtx > 10~ 6 are 
considered. 

inclusion of the charge conjugate mode is implied. 
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The primary background to the reconstructed decay channel arises from light quark continuum 
events e + e~ — > qq ( q = u,d,s,c ), which are characterized by collimation of final state particles 
with respect to the quark direction, in contrast to the more spherical BB events. Exploiting these 
different topologies we can increase the signal significance using topological variables computed from 
the center-of-mass (CM) momenta of all reconstructed charged and neutral particles in the event. 
For each event we linearly combine the sphericity, the angle between the B thrust axis and detector 
longitudinal axis, and the zeroth and second order Legendre moments into a Fisher discriminant (J-) 
[12] ■ whose coefficients are chosen to optimize the separation of signal and continuum background 
Monte Carlo samples. After optimization of the selection with respect to the simultaneous variation 
of all the selection criteria, we obtain that the Fisher requirement retains 71.2% of the candidates 
from the signal Monte Carlo sample and only 6.4% from the continuum background Monte Carlo. 

To further reduce combinatoric background we take advantage of the long mean lifetime of A 
particles and require that the separation of the A and B° vertices divided by its measurement error, 
computed on a per candidate basis by the fit procedure, exceeds 35. This criterion was optimized 
on Monte Carlo events and is effective in rejecting 42% of combinatorial background that survived 
all other cuts, while retaining 90% of signal candidates. The only sizable B background is from the 
process B° — > A c (— > Ak~) p, and we reject this with a veto on candidates with an invariant mass 
m(Air), within 20MeV/c 2 , approximately 5 standard deviations, of the nominal A c mass |10| . 

The kinematic constraints on B mesons produced at the T(4S) allow further background dis- 
crimination using the variables tties an d AE. We define m^s = v/(f + Pi • Pb) 2 /Ef — pb 2 where 
(Ei,pi) is the four momentum of the initial e + e~ system and ps the momentum of the reconstructed 
B candidate, both measured in the laboratory frame, and s is the square of the total available en- 
ergy in the T(4S) center of mass frame. We have AE = E B — where E* B is the B energy in the 
T(45) center of mass frame. Candidates satisfying \AE\ < 200 MeV and 5.2 < toes < 5.29GeV/c 2 
are used in the maximum likelihood fitting process. 

4 BRANCHING FRACTION 

We perform the measurement using a maximum-likelihood fit on ities-AE observables of recon- 
structed B candidates. The s "Plot technique ^3] , is then used to determine the m(Ap) distribution 
of reconstructed candidates and, once the correction for the nonuniform reconstruction efficiency is 
applied, measure the m(/lp)-dependent differential rate together with the total branching fraction. 

We consider as signal events only reconstructed B meson candidates whose daughters are cor- 
rectly assigned in the decay chain. By self-cross- feed, we refer to candidates reconstructed as signal 
events in which one or more particles are not correctly assigned in the decay chain. Examples of 
such misreconstruction include events in which a proton from the other B meson are associated to 
the signal B, and events where the protons from the signal B and A decays are interchanged. We 
define the total PDF in the AE-m^s plane as the sum of signal, self-cross-feed, and background 
components: 

1 N 
£ = — e-^+^B+iVs/^) -q [NsVs ^ + NBVBa + Nsfsc{Vsc{a] (i) 

where the product is over the N fitted events with iVg and N-q representing the number of signal 
and background candidates and / sc f representing the self-cross-feed fraction. The three V functions 
are taken as products of 1-dimensional AE and toes PDF's. We are justified in this simplification 
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by the small correlation between these two variables in our Monte Carlo sample, measured as -7.4% 
for signal, and -0.5% for background. The rngs PDF is taken as a double Gaussian for the signal 
and a threshold function |13j for the background. The AE PDF is taken as a double Gaussian for 
the signal and a first order polynomial for the background. Finally, the self-cross-feed contribution 
shows a peaking component that is modeled as the product of a double Gaussian in AE and a single 
Gaussian in ttt-es- The self cross- feed fraction / sc f = 0.59%, and the other parameters that enter 
the definition of its contribution to the PDF have been determined from a Monte Carlo sample of 
signal events. 

We vary the means of the narrow AE and mss signal Gaussians, the coefficient in the exponen- 
tial of the Argus function, the linear coefficient of the AE background distribution, and the event 
yields iVs and Nb- The means of the wide Gaussians are determined by applying Monte Carlo 
determined offsets to the mean of the narrow ones, such that only an overall shift of the fixed PDF 
shape is allowed. 

Once the maximum likelihood fit provides the best estimates of the PDF parameters, we use 
the s V\ot technique Jl] to reconstruct the efficiency-corrected m(Ap) distribution and measure 
the branching fraction. The PDF is used to compute the s-weight for the n-th component of event 
e as: 

s y n (y e) - E n LiNkVk{ye) W 

where the indices n, j and k run over the n c = 2 signal and background components whose dis- 
tributions, as functions of y e = (mEs^,AB e ), are identified with the Vj(y e ) symbol. V nj - is the 
covariance matrix of the event yields as measured from the fit of the PDF to the data sample. An 
important property of the s Vlot is that the sum of s-weights for the signal or background compo- 
nent equals the corresponding number of fitted signal or background candidates. Thus the gVlot is 
a good estimator of the m{Ap) distribution, and preserves the total signal yield, as determined by 
the maximum likelihood fit. To retrieve the efficiency-corrected number of events in given m{Ap) 
bin J we use the s-weight sum: 

Nj = ^^LLL. (3) 



sVn(y e 
seJ e{x e ) 



where e(x e ) is the per-event overall efficiency. The reconstruction efficiency depends on the position 
of the candidate on the square Dalitz plane x e = {m^ p , cos(#h))> $h being the helicity angle of the 
pion in the Ap rest frame, and has been measured on a 10 x 10 grid over the square Dalitz plane, 
using fully reconstructed signal Monte Carlo events. The error a [Nj] in Nj is given by: 

G , [Nj] = E ( (4) 



eeJ 



An estimate of the efficiency-corrected number of events in the sample is given by the sum of the 
efficiency-corrected s-weights or 

J 

and the total branching fractions is obtained from 

N 

B(B -» Apir) = — — r. (6) 
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Table 1: Systematic errors on the BF measurement. 





source 


error 


Overall 


BB counting 


1.1% 




B°B°/BB fraction 


1.4% 




Tracking efficiency 


3.9% 




PID efficiency 


1.4% 




MC statistics 


2.0% 




A — > pir branching fraction 


0.8% 


Event Selection 


Event shape cut efficiency 


2.4% 




Fit probability cut efficiency 


5.0% 




A flight length cut efficiency 


2.8% 




A mass cut efficiency 


2.4% 




A c veto cut 


0.5% 


Fit Procedure 


Likelihood parameters 


3.9% 




AE resolution 


1.7% 




Self cross-feed fraction 


0.8% 




s Vlot bias correction 


0.6% 


Total 




9.4% 



Using fully reconstructed signal Monte Carlo events, we have checked that this procedure pro- 
vides a measurement of the m(Ap) distribution and total branching fraction with negligible biases 
and accurate errors. 

5 SYSTEMATIC ERRORS 

Systematic errors are listed in Table I and classified as overall uncertainties, uncertainties associ- 
ated with event selection, and uncertainties associated with fitting the signal event distribution. 
Tracking efficiency uncertainty dominates the first category with a contribution of 3.9%. Particle 
identification systematic errors were evaluated by studying data versus Monte Carlo agreement of 
identification efficiency on protons from a pure sample of A — > pir decays and pions from K$ — > -kit 
decays. The finite signal Monte Carlo sample available to measure the reconstruction efficiency 
over the Dalitz plane, results in an additional 2.0% contribution to the systematic error. The 
uncertainty in the determination of the number of BB pairs in the data sample accounts for a 
1.1% systematic, while the assumption of a 50% ratio of B°B° to BB at the T(4S) gives an ad- 
ditional 1.4% contribution, computed as the difference with respect to the current measured value 
49.3 ±0.8% P33. 

Event selection systematic errors associated with the determination of the efficiencies of the 
Fisher-discriminant event shape cut and the vertex fit-probability cut, have been evaluated com- 
paring data and Monte Carlo selection efficiencies of a sample of B° — > J/tpK® candidates. In 
addition, we use an inclusive sample of A — > pir candidates to estimate systematic errors associated 
with the determination of efficiencies of flight length-significance cut and Tl-mass requirement. 

The application of the requirement on the reconstructed m(Air) invariant mass to veto B° — > A c p 
background has two associated systematic effects. The first causes an approximate 0.2% increase in 



11 



5.2 5.21 5.22 5.23 5.24 5.25 5.26 5.27 5.28 5.29 

(GeV/c 2 ) 



-0.2 -0.15 -0.1 -0.05 -0 0.05 0.1 0.15 0.2 

AE (GeV) 



Figure 1: Left plot: toes distribution of candidates with \AE\ < 27MeV. Right plot: AE distri- 
bution of candidates with toes > 5.274 GeV/ c 2 . Superimposed are projections of the 2-dimensional 
fit PDF onto the respective axes. 

the branching fraction due to the residual A c component that survives the cut. The second causes 
an approximate 0.5% reduction of the branching fraction due to the reduced Dalitz-plot space. We 
take the larger of the two as the systematic error associated with the A c veto cut. 

We vary parameters that are kept fixed in the likelihood fit by their statistical errors, as mea- 
sured on the signal MC sample fit, and measure the variation of the gVlot fitted result. The 
changes associated to the parameters that enter the definition of the signal PDF are conservatively 
considered as fully correlated and added linearly to give a signal PDF systematic error of 3.2%, 
where the uncertainty on signal toes fixed parameters accounts for a 1.9% contribution and the 
uncertainty on signal AE fixed parameters for a 1.3% contribution. The same procedure is applied 
to the parameters that enter the background PDF definition, with errors determined on luminosity- 
weighted background MC samples, giving an additional 2.2% systematic error. Finally, we combine 
in quadrature the two errors and obtain a 3.9% systematic error associated with uncertainties on 
the shape of signal and background PDF models. The comparison of B -> J/%I)K° S data and Monte 
Carlo samples reveals that the width of the AE Gaussian in the signal PDF can be underesti- 
mated in the Monte Carlo by up to 5%, and this translates to an additional 1.7% systematic error 
associated with the uncertainty in the AE resolution. 

We estimate possible biases associated with the determination of yields with the s V\ot tech- 
nique, using a collection of Monte Carlo experiments in which signal candidates, generated and 
reconstructed with a complete detector simulation, have been mixed with background candidates, 
choosing numbers of signal and background candidates similar to those expected on the data sam- 
ple under study. Biases have been found within the statistical error in their measurement, and we 
estimate a 0.6% systematic uncertainty associated with the s V\ot fitting technique. 

6 RESULTS 

We select a total of 4261 candidates in the region \AE\ < 200 MeV, to E s > 5.2 GeV/ c 2 , \m(Air) - 
m(A c )\ > 20 MeV/ c 2 in the 210.3 fb _1 data sample considered. Table II contains the fitted values of 
the 2-dimensional m^-AE PDF parameters, while Fig. ^ shows projections of the 2-dimensional 
PDF on the toes an d AE axes. Figure [2] shows the efficiency-corrected signal s V\ot distribution 
of candidates as a function of the m(Ap) coordinate; this reveals a near-threshold enhancement 
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Figure 2: s "Plot of the m(ylp) event distribution with efficiency corrections applied. 

Table 2: Likelihood fit result. N$ and Nb are the number of fitted signal and background candidates 
respectively. \i (AE) is the mean value for the narrow Gaussian of the AE signal PDF component, 
while ci (AE) is the slope of the linear AE background PDF. fx (tties) is the mean value for 
the Gaussian of the toes signal PDF, and CArgus ("t-es) is the coefficient at the exponent of the 
background toes Argus function as given in Reported errors are statistical only. 



Parameter 


Value 


N s 


7 q 7 +12.U 

'J- 1 -11.2 


iV B 


4187 ± 66 


H(AE) 


-1.71 ±3.10 MeV 


d (AE) 


-3.71 ±0.25 


V( m Es) 


5.2808 ± 0.0004 GeV/c 2 


CArgus (WES) 


-15.1 ± 1.7 



similar to that observed in other baryonic B decays. Summing the efficiency-corrected s"Plot bins, 
we obtain a yield of 488 ± 79 signal events, where the error is statistical. Using Equation El we 
measure the branching fraction: 

B(B° -> Apir-) = [3.30 ± 0.53 (stat.) ± 0.31 (syst.)] x IO" 6 . 

This measurement, which is compatible with a previous measurement by the Belle collaboration^, 
confirms the peaking of the baryon-antibaryon mass spectrum toward threshold, a feature that plays 
a key role in the explanation of the higher branching fraction of three-body baryonic B decays with 
respect to two body ones. 
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